Mon. Not. R. Astron. Soc. OOP. [TlfTil (20111 Printed January 12, 2013 (MN 1^1^ style file v2.2) 



Steady-state evolution of debris disks around solar-type 
stars 

N. Kainsi'2 ^ M.C. Wyatt^ J.S. Greaves^ 

^SUPA, School of Physics and Astronomy, University of St. Andrews, North Haugh, St Andrews, KY16 9SS, United Kingdom 
^European Southern Observatory, Karl-Schwarzschild Strafie 2, 85748 Garching bei Miinchen, Germany 
^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge GB3 OHA, United Kingdom 



Accepted ... Received ... ; in original form 



ABSTRACT 

We present an analysis of debris disk data aroun d Solar-type s tars (s pectral types 
F0-K5) using the steady-state analytic al model of Wvatt et al.l (|2007a[ ). Models are 



fitted to published data from the FEPS (jMever et al.ll2006b[ ) project and various GTO 
programs obtained with MIPS on the Spitzer Space Telescope at 24^m and 70^m. and 
compared to a previously published analysis of debris disks around A stars using the 
same evolutionary model. We find that the model reproduces most features found in 
the data sets, noting that the model disk parameters for solar- type stars are different 
to those of A stars. Although this could mean that disks around Solar- type stars 
have different properties from their counterparts around earlier-type stars, it is also 
possible that the properties of disks around stars of different spectral types appear 
more different than they are because the blackbody disk radius underestimates the true 
disk radius by a factor Xr which varies with spectral type. We use results from realistic 
grain modelling to quantify this effect for solar-type stars and for A stars. Our results 
imply that planetesimals around solar-type stars are on average larger than around A 
stars by a factor of a few but that the mass of the disks are lower for disks around 
FGK stars, as expected. We also suggest that discrepancies between the evolutionary 
timescales of 24^m statistics predicted by our model and that observed in previous 
surveys could be explained by the presence of two-component disks in the samples 
of those surveys, or by transient events being responsible for the 24//m emission of 
cold disks beyond a few Myr. Further study of the prevalence of two component disks, 
and of constraints on Xr, and increasing the size of the sample of detected disks, are 
important for making progress on interpreting the evolution of disks around solar-type 
stars. 
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1 INTRODUCTION 

Recent surveys looking for excess emission around solar- 
type (F, G and K spectral types, hereafter FGK) stars 
have found that ove r 15% of thes e objec ts are surrounded 



by d eb ris disks jBrvden et ah l2006l: Beichman et al.l 



2006a; Moro-Martm et al 
Hillenbrand ct al 



2008 



2007bl: iTriUing et all l2008l : 



Greaves et al.l l2009h . Although 



these disks are believed to be analogues of the Solar System's 
asteroid and Kuiper belts, the infrared luminosity of most 
detected disks around solar-type stars is usually around 100 
times higher than that of the Solar System's debris struc- 
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tures jGreaves fc Wyattll2010ll . This suggests that the Solar 
System may have unusually sparse debris, or that the de- 
tected disks around other stars are anomalously luminous, 
perhaps as a result of recent collisions in their planetesimal 
belts, producing transient excess dust. This may be caused 
by a process similar to the Late Heavy Bombardment (LHB) 
perio d which occured ~ 700 Myr after the formation of the 
Sun (|Booth et al.ll2009l '). 

Among the sample of FGK stars with detected disks, 
several are known to harbour extrasolar planets. Inter- 
estingly, searches for asteroid and Kuiper belts analogues 
around some of the known extrasolar planetary systems or- 
biting solar-type stars, such as r Boo, tt Men (HD 39091), 
V And, 55 Cnc and 51 Peg did not find detectable lev- 



els of dust (|Brvden et all l2009l : iBeichman eFall l2006lj) . 
Larger-scale searches for a correlation betw een the pres- 
ence of a planet and that of a debris disc jKospal et aP 
I2OO9I : iMoro-Martin et al.l l2007al : iBrvden et all |2009| ) also 
concluded that the rates of excess emission were not sig- 
nificantly different when comparing stars with and without 
planets. These results could indicate that dust luminosity 
has fallen to near the luminosity of the Kuiper and aster- 
oid belts of the Solar System and that these systems are 
past their own LHB equivalent. Such a system would not 
only not feature the transient excess dust associated with 
such an event, but the fact that cleari ng occurs as part o f 
the process that leads to the LHB (e.g. iMever et al.ll2006al ) 
means that dust e xcess emission would have to be faint. 
iBooth et ai] (|2009l ). however, suggest that such events are 
rare, and an alternative explanation could be that these sys- 
tems might also have been born with lower-mass planetesi- 
mal belts. Studying extrasolar systems offers an important 
perspective for understanding the history of the Solar Sys- 
tem. In particular, it is essential to gain an understanding of 
how planetesimal belts evolve around their host stars, and 
which mechanisms affect their structure and their luminosity 
in order to be able to compare observed extrasolar systems 
to t he Kuiper belt. 

iKenvon fc Bromle^ (|201(]| ') published calculations of the 
formation of icy planets and debris disks at 30-150 AU 
around 1-3M@ stars and found observations of A- and G- 
type stars to favour models w ith small {'--^ 10km) planetesi- 
mals. ICarpenter et al.l (|2009l ) analysed the evolution of cir- 
cumstellar dust around solar-type sta rs using the model o f 
iDominik fc DecinI (|2003l '). extended bv lWvatt et~all (|2007al ). 
This simple analytical model of steady-state debris disk evo- 
lution was developed by considering the coUisional grinding 
down of planetesima l belts. Here we also use the model of 
IWvatt et all |2007i) to carry out an analysis of the statis- 
tics of excess emission around FGK stars. The basic fea- 
tures of the model are recalled in Sec. [2] In Sec. [3l we 
outline our modelling procedure. In Sec. 31 we briefly de- 
scribe the 24 and 7 0/^m data wc us e, taken from the stud- 



ies of iTrilling et al. (2008 1. Hillenbrand et al.l (|2008l ') and 
IBeichman et all (|2006al ). In Sec. [3 we fit this data with our 
model and discuss our results in Sec. [SI comparing them to 
the results of the analogous stu dy that was carried out for 
A stars by IWvatt et ail (|2007bt l. We also discuss the pos- 
sible reasons for differences between the evolution of debris 
disks around FGK stars and that of their earlier-type coun- 
terparts, and consider individual systems that are not well 
fitted by our models. We conclude by looking at the implica- 
tions of our models for the properties of disks around FGK 
stars compared to those around A stars. 



2 DEBRIS DISK MODEL 

In this section we re call the main features of the analytical 
model developed by IWvatt et al.l (|2007al ') and reprised by 
IWvatt et al.l (|2007bl ) . as well as the assumptions made in 
applying this model to populations of optically thin debris 
disks. We start with a planetesimal belt characterised by a 
size distribution 



where D is the diameter of the planetesimals; this is the 
distribution expected for a planetesimal belt in coUisional 
equilibrium (|Dohnanvilll969l ). This distribution is assumed 
to be valid from the largest planetesimals, of diameter Dc 
down to the blowout diameter Dbi, below which particles are 
blown away by radiation pressure. 

For a single-radius planetesimal belt at radius r, with 
a width dr, the fractional luminosity of the dust emission, 
/ = Lm/Lt can be expressed in terms of the total cross- 
sectional area crtot as 



f 



O"tot 



(2) 



Given the size distribution in Eq. Ill, crtot is propor- 
tion al to Mtot through a constant that depends on Dbi and 
Dc (|Wvatt et al.ll2007al ). 

Assuming the dust particles act as blackbody emitters, 
the blowout diameter is given by 



0.8 



L, 2700 



(3) 



where L, and M* are in solar units, and p is the density of 
the dust particles in kg m~^. The temperature T of the dust 
at radius r is given by 



T = 278.3 L°-^\'°-^ . 



(4) 



For blackbody emitters, we also know that emission 
from the disk at wavelength A is 



F,.,d^sk = 2.35 X 10 "B„(A,r) crtot d~ 



(5) 



where d is the distance to the star in pc and is in Jy if 
the Planck function is in Jy sr^^. 

Considering a coUisional cascade in which the popula- 
tion within a given size range is being destroyed in collisions 
with other members of the cascade, but is being replenished 
by fragmentation of larger objects, the collision of the larger 
objects in the cascade solely determine the long-term evo- 
lution of the population. Hence the long-term evolution de- 
pends of the coUisional lifetime tc of the planetesimals in the 
disk with size D^. The expression for tc can be expressed as 



3.8pr^-^{dr/r)Dc 



(6) 



n{D) <x D~ 



(1) 



where tc is in Myr, dr/r is the fractional width of the disk, 
Dc is in km, Mtot is in units of A/0, e and / are the eccen- 
tricities and inclinations of the planetesima. ls ' orbit s, and 
the factor G(Xc) is defined in IWvatt et al.1 l|2007ah : note 
that this simplified expression is only valid when e — I 
(which we will assume throughout this paper). Xc is de- 
fined as Xc = Dcc/Dc, where Dec is the diameter of the 
smallest planetesimal that has sufficient energy to destroy a 
planetesimal of size Dc- This factor can be calculated from 
the value of the dispersal threshold , which is defined as 
the specific incident energy required to catastrop hically de- 
stroy a particle. It follows (|Wvatt fc Dentll2002l ') that with 
Q*D given in Jkg"\ 
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1.3 X 10" 



iD 



rM: 



2.25 e2 



1/3 



(7) 



Ignoring non-coUisional processes, the time-dependence 
of the disk mass can be worked out by solving the differential 
equation dMtot/dt = —Mtot/tc, yielding 

Mtot(i) = Mtot(0)/[1 +fAc(0)]. (8) 

As noted by I Wvatt" et all (|2007a!), the mass of the disk 
at t ^ does not depend on the initial disk mass, since tc{0) 
depends on Mtot(O). As a result, at a given age, there is a 
maximum mass that can remain after coUisional evolution, 
and therefore a maxim um infrared luminosity /max, given 
by (|Wvatt et al.ll2007ah 



/"max — 



10" 



\dr/r) 



for Dbi in ^m. 

For observations limited by a calibration limit expressed 
as a flux ratio -Rdct(A) = -F^,disk/-Piy,phot, the corresponding 
detection limit in terms of fractional dust luminosity, /dot, 
is given by 

/dot = 6xl0''i?dctr-'L.T-''B,.(A,r.)/S,.(A,278.3LrV"°'; 

(10) 

This is the threshold value we use to determine which 
model radii are detectable at given wavelengths in the rest of 
this paper. Finally, the limit at which Poynting-Robertson 
(PR) drag becomes important, i.e. when radiation and gravi- 
tational timescales become s imilar for the smallest particles, 
is given l|Wvatt et aLlbOGTal ) by 



/pR = 50 X 10"" {dr/r)^hh/r . 



(11) 



That is, if a star system has a flux / < /pr, then it is 
likely that it will be affected by PR drag, which could mean 
that in such a disk the dust component of the disk becomes 
spatially separated from the planetesimal belt. 



3 MODELLING PROCEDURE 



As done bv I Wvatt etlo] (|2007bh for A stars, we apply this 
model to debris disk populations assuming that all stars 
have a planetesimal belt which is undergoing collisional pro- 
cesses as described by the above equations. This means that 
the initial state of the belts is completely determined by 
the parameters r, dr, p, Dbi, Mtot(0),L, and M, (of which 
the latter two determine Dbi through Eq. |3]), while their 
evolution also depends on Q*jj,e,Dc and I. We also set 
p = 2700kg m-^,dr = r/2 and e = / = 0.05. Values of 
Q*]j and Dc were then constrained by fits to observational 
data of debris disks. 

Values of A/tot (0), r, L, and M* were drawn from dis- 
tributions of properties of stars and debris disks. We used a 
power law distribution for disk radii, N{r) oc r'', and treated 
7 and the limits of the radius distribution, rmin and rmax, 
as free parameters. 



mode 

Trilling et ol. (2008) 
Hillenbrond et ol. (2008) 
Beictimon et ol. (2006) 



Figure 1. The distribution of spectral types for the model pop- 
ulation, obtained by combining the observed distributions from 
the different surveys. The individual distributions are also plotted 
in red dotted line (Trilling et al. 2008), blue dashed line (Hillen- 
brand et al. 2008) and green dash-dot line (Beichman et al. 2006). 



The distribution of disks at both wavelengths is affected 
by strong observational bias, especially for solar-type stars. 
For these stars, disks with radii larger than ~ 50 AU can 
only be detected at 24/im if they have a large fractional ex- 
cess (see e.g. Fig. |4] and [8]), meaning that there is a bias 
towards the detection of disks with radii lower than ~ 50 
AU, and this must be accounted for by a satisfactory model. 
Therefore we constrained 7 by comparing the fit of the ob- 
served distribution to a fit to the subsample of the model 
population which could be detected at 24 and 70/^m accord- 
ing to Eq. (|10l) . rather than a fit to the whole model popu- 
lation. This comparison also yielded constraints on T'niin and 
''max, although these parameters are also sensitive to the 
fraction of disks detected in different age bins. 

For the distrib ution of initial dust disk masses M tot (0) , 
we use the resuhs of lAndrews fc Wifliamj (|2005l . l2007l 'l. who 
derived a lognormal distribution of dust masses centred on a 
value Mmid, with a standard deviation of 0.8 dex, from sub- 
millimetre observations of protoplanetary disks in Taurus- 
Auriga, i.e. for ~ solar-mass stars. We use their value for the 
distribution's dispersion but fit Mmid as a free parameter. 
We choose this approach because the submillimetre data of 
these studies do not detect all sizes of dust grains. 

A spectral type was drawn on the range FO - K5, chosen 
to correspond to the range observed by the various Spitzer 
programmes jTrilling et all I2OO8I : iHillenbrand et all (20081 : 
iBeichman et al.ll2006ai r The distributions of spectral types 
observed in these programmes were used as distributions for 
generating our model population; the resulting model distri- 
bution is plotted on Fig. [T] We then determined a value of 
Dbi for each spectral type, using corresponding values of L* 
and M* and Eq. ((3]). We also assigned values to other pa- 
rameters using the relevant distributions, determining the 
initial conditions of each disk. It should be noted that this 
approach assumes that the mass and radius of the disk are 
independent both of each other and of the properties of the 
star around which they are located. 

The systems were then evolved using Equations [6] and 
[S] and values were drawn at random for the age and dis- 
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Trilling et ol. (2008) 
Hillenbrond et ol. (2008) 
Beichmon et ol. (2006) 



Age [Gyr] 

Figure 2. The distribution of stellar ages for the model popula- 
tion, obtained by combining the observed distributions from the 
different surveys. 
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Trilling et ol. (2008) 
Hillenbrond et ol. (2008) 
Beictimon et ol. (2006) 
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distance [pc] 

Figure 3. The distribution of distances to the stars for the model 
population, obtained by combining the observed distributions 
from the different surveys. 



tance of each system, in order to determine its current "ob- 
served" properties. Distributions for these properties were 
also chosen to reflect the the data sample we used (Fig. [2]- 
[3| . The range of values for age and distance were between 
and 10 Gyr, and between and 150 pc respectively, with 
most stars between and 60 pc (chosen to cover the same 
ranges as the observational data used in this paper). We 
do not take into account any correlation between age and 
distance, but although we assign a distance to each system, 
its value has no impact on the flux statistics on which the 
analysis presented in this paper is based, because all ob- 
servations are assumed to be calibration-limited; we do not 
treat the distribution of upper limits (see Sec. U). This left 



QD,-Dc,7,Mmid, 



mm, ' max 



as free parameters of the model. 



Finally, a good-fit model was only retained if the parameter 
values were realistic, i.e. consistent with ra nges of values pre- 
dicte d by models of catastrop hic collisions (|Benz fc Asphaud 
Il999l) and planet formation (jKenvon fc BromlevlbOOa V 



4 DATA 

A list of all stars with 70^m excess emission only in the 
sample we use in this paper is given in Table [TJ while a list 
of sources that show excess at both 24 and 70/im is given 
in Table (2] stellar properties for these are published values 
resulting from Kurucz model fits, or in cases for which these 
were no t available, were computed using Schmidt-Kaler re- 
lations l|Aller et al.lll982l ). Using the 24 - 70/im colour tem- 
perature of the dust, and assuming that the dust acts as a 
blackbody, we also calculated a disk radius for the disks in 
Table [2] with Eq. (Bll . We do not include the debris disks 
published recently bv lKoerner et al.1 (|2010l ) because they do 
not report stellar ages or non-detections in their data. This 
sample would only add one disk detected at both 24 and 
70/im, and therefore would not infiuence our results signifi- 
cantly. 

We made a cut in stellar age, excluding objects with 
an age below 30 Myr in order to avoid protoplanetary disks 
affecting our statistics. We also excluded disks for which 1-cr 
error bars were above a threshold at — 5 mjy. This thresh- 
old was chosen empirically to avoid having noisy fiux upper 
limits being wrongly counted as large exce sses. These cut s 
affected mostly the FEPS data ( Hillenbr and et al.l [ioosi '). 
with 6 disks detected at both wavelengths being removed 
from their published sample. The reason for this is the dif- 
ference in the observational approache s of the SIMTPF and 
FE PS observati ons: the observations of lTriUing et al.1 (|2008l ) 
and lBeichman et al. (2006a) at 70/j m are c alibration-limited 
while those of iHillenbrand et al.l (|2008l ) are sensitivity- 
limited; at 24/im, all observation are calibration-limited (see 
Table [21 . The consequence of this is that most of the FEPS 
disks are very bright disks, since these are the only ones for 
which good photometry could be obtained. 

The resulting sample consists of 46 disks detected at 
70/im only and 17 detected at both 24 and 70/^m. Below is 
a short description of each subsample we used. 



4.1 Trilling et al. (2008) sample 

iTriUing et akl (|2008l ) collected a sample of 193 F, G and K 
stars observed with the Multiband Imaging Photometer for 
Spitzer (MIPS). To gether with data already published by 
iBrvden et al.l (120061 ). the whole sample gives a view of de- 
bris disks around stars with masses and ages similar to that 
of the Sun, covering a range of spectral types between FO and 
K5. Their "solar type" sample was selected using criteria on 
the spectral type (F5-K5) and luminosity class (IV or V), 
as well as setting a minimum photospheric 70/im fiux and 
signal-to-noise (S/N) ratio. From the combined FGK sam- 
ple, 27 show significant excess at 70/im, meaning an excess 
with X70 > 3, where X70 is defined as 



X70 = 



F7 - P70 

O"70 



(12) 



where F70 is the total 70/im fiux measured, P70 is the 
predicted photospheric flux at 70/im, and a"7o is the error 
bar associated with the fiux measurement. On top of these 
27 objects, 3 m ore are classifi e d as e xcess sources for rea- 
sons detailed in iTrilling et al] (|2008l ). but one is removed 
due to the stellar age not being determined, bringing the 
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total number of excess objects in the sample to 29. F inally, 
a re-reduction of data for HD 101259 jSu et al.ll2oTol ') found 
that this system did not in fact have significant excess, so 
we removed this object from the sample. This left 28 objects 
with excess emssion, including 6 systems for which signifi- 
cant excess flux is detected at both 24 and YOfim. 



4.2 Beichman et al. (2006) sample 

The sample of iBeichman et all (j2006al ) includes objects ob- 
served in the frame of other projects, including radial ve- 
locity search teams, coronagraphy and interferometry mis- 
sions. As a result, it includes some low-mass close stars. 
They selected stars within 25 pc (with a few exceptions 
for the earlier-type stars) and excluded targets with binary 
companions within 100 AU on the grounds that binarity 
might prevent the formation or long-term stability of plane- 
tary systems. Their selected sample was then observed with 
MIPS at 24 and 70^m and four stars were observed further 
with the Infrared Array Camera (IRAC) in order to help 
determine their photospheric fiux. 

Based on the c riterion expressed by Eq. (|12p . 
IBeichman et al.l (|2006al ) identify 12 stars out of the 88 in 
their sample to have a 70/im excess at the 3 a level or better. 
Amongst these, 5 also display significant excess at 24^m. 



— <> — < _ 
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Figure 4. Histogram for the distribution of disk radii of the 
combined observational samples (diamonds) for disks detected at 
both 24 and 70 fim, shown with error bars calculated for small 
samples by Gehrels et al. (1986). Also plotted are the distributions 
for the entire model population (red, dashed line) and the model 
population that could be detected at 24 and 70/im (blue, dotted 
line). 



4.3 Hillenbrand et al. (2008) sample 

iHillenbrand et ahl (|2008l ) presented data obtained using 
MIPS, IRS and IRAC, with most of their systems ob- 
served at 3.6, 4.5, 8.0, 13, 24, 33, 70 and 160/im, as part of 
the FEPS {Formation and Evolution of Planetary Systems) 
program. They report 25 systems with excess flux at 70/im 
and S/N.^y^j^ > 3. Amongst this sample, 13 also have excess 
emission at 24/im. Since photospheric sensitivity at 70/xm 
could not be achieved without very long integration times 
except for the closest stars in the sample, the sensitivity 
of their observations was determined by a target detection 
threshold of dust emission, which they expressed relative to 
dust emission in a young Solar System model. For most tar- 
gets, the survey was sensitive to 5-10 times the dust emission 

predicted by the model. 

The sample oflHillenbrand et all (120081') is different from 
those of lTrilling et all (|2008h and IBeichman et all (|2006al '). 
as the former sample is age-selected, leading to an even 
distribution in logarithmic age bins, while the latter are 
volume-limited, leading to a linear age distribution. This 
difference is taken into account in the analysis that follows 
when fitting power laws to time-dependent parameters, and 
the data cuts we made were chosen to minimise the effect of 
poorly constrained photospheric fiuxes in the FEPS sample 
on our modelling of the statistics. 



5 RESULTS 

5.1 Fit to the 70/im statistics 

Following the procedure described in Sec. O we found 
best-fit values of Dc = 450 km and = 3700 J kg' \ 
These are consistent wi th models of catastrophic collisions 
jBenz fc Asphauelll999l ). We also found best-fit parameters 



for the radius distribution of rmin = 1 AU, rmax = 160AU 
and 7 = —0.60 ±0.35, simila r to the value of 7 ~ —0.8 =fc 0.3 
found in the A stars study bv lWvatt etall lj2007bt ). The val- 
ues of a were obtained using error bars on measured fluxes 
reported in the literature. A histogram of the distribution of 
disk radii for the model and observed data is shown on Fig. ID 
In order to reflect our small sample sizes, we use er ror bars 
for sm all-number Poisson statistics as calculated bv lGehre3 
()l986l ) . Our radius distribution exten ds to lower radius than 
the one found bv lLohne et al.l (|2008h . who used a value of 
rmin ~ 20 AU. However they limited their data sample to G 
stars, and their analysis used radius values calculated by as- 
suming realistic grain emission rather than blackbody emis- 
sion. Furthermore, we argue in Sec. I6.2l that the true radius 
distribution should be larger than the one we derive. Fi- 
nally, we find a best-fit value for the median of the initial 
dust mass distribution of Mmid = 4M0, which is con sistent 
with r esults for the disk mass-stellar mass relation of iNattal 
(|2004h . who find an approximate range of total disk masses 
around solar-mass stars of —2.5 < log(Mdisk/Af«) < —0.5. 
This val ue is lower than that of Mmid ~ WM^ found for A 
stars bv lWvatt et al.l (|2007bl ). 

The models found with these parameters are plotted on 
Fig. O which show that the statistics are repro duced by the 
model convincingly, as they were as for A stars (|Wvatt et al.l 
l2007bl ): this is shown on the right-hand panel of Fig. [S] The 
fiux ratio shows very slow evolution past the earliest age 
bins, which is what is seen in the data as well. We use two 
excess categories, small (i?70 < 15) and large {R70 > 15). 
The threshold value i?70 = 15 was chosen empirically to 
avoid upper limits contaminating the sample of stars with 
large excess emission. The large excess fraction makes up 
~ 20% of excesses at early ages, falling to a few percents 
within ~ 3 Gyr. 
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Table 1. List of sources with 70 fim excess only. Flux ratios and limits are given as total flux divided by photospheric flux. i?24 is 
deflned as _R24 = F24/F24., and similarly, _R70 = F70/F70,. C24 and C70 are the 3 - calibration limits at 24 and 70/xm, deflned as 
Ca = l + 3a.,A/i^.,A. 



Star name 


Sp. type 


d 


iage 


R2A 


C24 


i?70 


C70 






(pc) 


(Gyr) 




{3<x) 




{3<x) 


Trilling ct al. 


HD 1581 


F9 V 


8.6 


3.02 


0.98 


1.10 


1.4 


1.3 


HD 3296 


F5 


47 


2.5 


1.02 


1.10 


4.8 


2.8 


HD 17925 


Kl V 


10 


0.19 


1.05 


1.10 


3.4 


1.7 


up* 1 QQQA 


-T O V 


99 


O.OO 




1 in 


1 7 


1.0 


HD 20807 


Gl V 


12 


7.88 


1.04 


1.10 


1.8 


1.5 


HD 22484 


F8 V 


17 


8.32 


1.02 


1.10 


1.9 


1.3 


HD 30495 


Gl V 


13 


1.32 


1.06 


1.10 


5.5 


1.6 


HD 33262 


F7 V 


12 


3.52 


1.03 


1.10 


1.7 


1.4 


HD 33636 


GO 


29 


3.24 


1.00 


1.10 


7.0 


2.2 


HD 50554 


F8 V 


31 


4.68 


1.00 


1.10 


8.4 


3.4 


HD 52265 


GO 


28 


6.03 


0.99 


1.10 


4.8 


2.9 


HD 57703 


F2 


44 


2.3 


1.03 


1.10 


9.3 


3.3 


HD 72905 


G1.5 V 


14 


0.42 


1.06 


1.10 


2.3 


1.5 


HD 75616 


F5 


36 


4.8 


1.03 


1.10 


10 


2.5 


HD 76151 


G3 V 


11 


1.84 


1.03 


1.10 


2.4 


1.6 


HD 82943 


GO 


27 


4.07 


1.02 


1.10 


17 


3.1 


HD 110897 


GO V 


17 


9.7 


0.98 


1.10 


4.3 


1.9 


HD 115617 


G5 V 


8.5 


6.31 


1.04 


1.10 


4.0 


1.5 


HD 117176 


G5 V 


18 


5.37 


0.98 


1.10 


1.8 


1.3 


HD 128311 


KO 


167 


0.39 


0.94 


1.10 


3.0 


2.3 


HD 206860 


GO V 


18 


5.00 


1.03 


1.10 


2.3 


1.5 


HD 212695 


F5 


51 


2.3 


1.00 


1.10 


9.5 


3.3 


Hillenbrand et al. 


HD 6963 


G7 V 


27 


1.00 


1.05 


1.13 


13 


8.5 


HD 8907 


F8 


34 


0.32 


1.05 


1.13 


46 


12 


HD 31392 


KO V 


26 


1.00 


1.02 


1.12 


20 


8.5 


HD 35850 


F7/8 V 


27 


0.03 


1.14 


1.14 


4.9 


3.9 


HD 38529 


G8 HI/IV 


42 


3.16 


0.96 


1.12 


4.3 


3.1 


HD 72905 


G1.5 


14 


0.10 


0.84 


1.10 


2.7 


2.1 


HD 122652 


F8 


37 


3.16 


1.08 


1.13 


24 


10 


HD 145229 


GO 


33 


1.00 


1.09 


1.14 


20 


9.1 


HD 150706 


G3 V 


27 


1.00 


1.05 


1.13 


8.7 


6.4 


HD 187897 


G5 


33 


1.00 


1.03 


1.12 


14 


7.5 


HD 201219 


G5 


36 


1.00 


1.07 


1.13 


18 


11 


HD 209253 


F6/7 V 


30 


0.10 


1.14 


1.14 


14 


6.8 


Beichman et al. 


HD 38858 


G4 V 


16 


4.57 


1.00 


1.30 


10 


3.0 


HD 48682 


GO V 


17 


3.31 


1.00 


1.59 


12 


2.0 


HD 90089 


F2 V 


21 


1.78 


1.00 


1.59 


2.3 


1.7 


HD 105211 


F2 V 


20 


2.53 


1.01 


1.19 


11 


2.4 


HD 139664 


F5 IV-V 


18 


0.15 


1.10 


1.24 


5.9 


1.6 


HD 158633 


KO V 


13 


4.27 


0.90 


1.53 


18 


2.0 


HD 219623 


F7 V 


20 


5.06 


1.04 


1.12 


3.0 


1.7 



5.2 Fit to the 24/im statistics 

Fig. m shows the model fit (with the parameter values given 
in the previous section) to the observed 24/im statistics. The 
statistics both in the model and the data samples show that 
24/im excess flux evolves on a timescale of ^ 1. — 1.5 Gyr, 
similar to what was found bv lLohne et al. | (|2008l). but slower 
tha n what is seen in the obs ervations of Isiegler et al.l (|2007l ) 
and ICarpenter et al.l (|2009l ) (see also Fig. |7]). 24/tm statis- 
tics are well fitted by the model and are particularly crucial 



in constraining the parameters which determine the evolu- 
tion timescale of the models (Qjj, e, Dc, Mmid), as they show 
clearer time evolution than the 70/im data; this is clearly 
visible on the left panel of Fig. |6] 

Fig. [7] shows the temporal evolution of the fraction of 
systems in the model population that have excess emission 
at 24/jm. Also p l otted for c omparison are the data presented 
bv lSiegler et al.l (|2007l ') and lCarpenter et all (|2009t) (t he lat- 
ter an extension of the results of Meyer et al.ll2008l). Ou r 
model is within \-2a of the results of Siegler et al.l (|2007f ) 
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Table 2. List of sources with 24 and 70/xm excesses. Flux ratios and limits are given as total flux divided by photospheric flux. R2A, 
R70, C24 and C70 as defined as in Table [T] Stellar properties are either published values found by fitting Kurucz models to the data 
when these were available, or calculated using Schmidt-Kaler relations for main sequence stars. 



Star name 


Sp. type 


d 


^age 




Af, 


/ 


r 


// fmax 


///pR 


R24: 


C24 


R70 


C70 






(pc) 


(Gyr) 


/Lq 


/Mq 


{/10-5) 


(AU) 








(3a) 




(3a) 


Trilling et al. 


HD 166 


KO V 


14 


5.0 


0.42 


0.79 


5.9 


9.1 


1.99 


8.0 


1.14 


1.10 


6.9 


1.8 


HD 3126 


F2 


42 


3.5 


2.9 


1.5 


13 


21.8 


1.86 


20 


1.16 


1.10 


27 


3.3 


HD 10647 


F9 V 


17 


6.3 


1.8 


1.1 


34 


21 


5.97 


59 


1.21 


1.10 


51 


2.1 


HD 69830 


KO V 


13 


4.7 


0.42 


0.79 


20 


1.0 


836 


9.0 


1.47 


1.10 


1.5 


1.5 


HD 105912 


F5 


50 


1.8 


3.2 


1.4 


7.9 


7.7 


5.71 


7.4 


1.52 


1.10 


11 


3.0 


HD 207129 


GO V 


16 


5.8 


1.5 


1.1 


12 


15.3 


3.42 


18 


1.17 


1.10 


16 


2.8 


Hillenbrand et al. 


HD 25457 


F7 V 


19 


0.10 


2.1 


1.3 


10 


17 


0.05 


14.8 


1.31 


1.16 


18 


5.0 


HD 37484 


F3 V 


60 


0.10 


3.5 


1.5 


32 


17 


0.23 


43.5 


1.43 


1.29 


46 


14 


HD 85301 


G5 


32 


1.0 


0.71 


0.92 


13 


8 


1.64 


14.8 


1.36 


1.17 


13 


8.4 


HD 202917 


G5 V 


46 


0.03 


0.66 


0.92 


25 


9 


0.08 


31.0 


1.63 


1.20 


28 


16 


HD 219498 


G5 


150 


0.32 


4.9 


0.92 


20 


31 


0.11 


46.4 


1.29 


1.15 


25 


14 


Beichman et al. 


HD 25998 


F7 V 


21 


0.6 


2.4 


1.2 


4.5 


13 


0.27 


5.8 


1.14 


1.12 


4.2 


2.2 


HD 40136 


Fl V 


15 


1.3 


4.3 


1.6 


1.9 


5.9 


2.24 


1.5 


1.13 


1.12 


1.6 


1.5 


HD 109085 


F2 V 


18 


1.3 


2.9 


1.5 


15 


5.9 


14 


11.8 


1.99 


1.12 


5.9 


1.6 


HD 199260 


F7 V 


21 


3.2 


2.4 


1.2 


3.3 


14 


0.96 


4.4 


1.11 


1.12 


3.5 


2.0 


HD 219482 


F7 V 


21 


6.1 


2.4 


1.2 


3.6 


18 


1.04 


5.5 


1.08 


1.12 


4.4 


1.7 
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Figure 5. Left: total 70/xm flux divided by the photospheric flux against age. The model population is shown with small dots, and the 
observations of Trilling, Hillenbrand and Beichman are shown as filled dots, triangles and squares respectively. Disks detected at TOfim 
with < 3(7 confidence are shown with open diamonds. Red symbols indicate disk emission detected at > 3cr at both 24/im and 70^m. The 
horizontal dashed line separates the populations with small and large excess fiux. Right: fractional populations of stars with different fiux 
ratios for different age bins (0.03 — 0.3 Gyr, 0.3 — 2 Gyr, 2 — 4 Gyr, 4 — 10 Gyr). Small (R70 = ^70,tot/i^70,phoi < 15) and large excesses 
(Rro > 15) are shown as crosses and squares respectiv ely. Model predictions are connected with lines and observed values l|Trilling et al.l 
I2OO8I : [Hillenbrand et al1l2008l : [Beichman et al.ll2006al ) arc plotted with small-sample Poisson statistics error bars (Gehrels et al., 1986) 
where appropriate. On both panels, the dotted vertical lines indicate the limits of the age bins. 



for early-age (< 100 Myr) systems, as well as those of 
[Carpenter et al.l ([20091 ') for systems with age > 100 Myr, 
but the temporal evolution we predict is significantly slower 
than that seen by those surveys. We discuss possible reasons 
for this in Section [6.31 

The combined fit to the 24 and 70^m statistics has a 
goodness-of-fit of x^/d.o.f. = 1.05. 



5.2.1 Discussion of individual sources 

From Fig. [6] it is clear that a number of systems are not well 
fitted by our models. We discuss shortly the systems labeled 
on the plot: HD 109085, HD 69830 and HD 10647. 

HD 109085 was found bv [Sheret et^ ([2004[ ') to have a 
radius of ~180 AU, compared to the blackbody radius of 6 
AU used in this analysis. This system was also reported to 
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Figure 6. Left: total 24/im flux divided by tlie pliotospheric flux against age. Symbols arc the same as for Fig. [5] A iiorizontal dasiied 
lines separates the populations with small and large excess flux. Right: fractional populations of stars with different flux ratios for 
different age bins (0.03 — 1 Gyr, 1 — 3 Gyr, 3 — 13 Gyr). Small (i?24 = ^24, phot < 1-25), and large excess es (R24 > 1.25) are 
shown as crosses, triangles and squares respecti vely. Model predictions arc connected with lines and observed values l l Trilling et al ] |2008l : 
iHillenbrand et al.|[20ol : iBeichman et al.ll2006a^ are plotted with small-sample Poisson statistics error bars (Gehrels et al., 1986) where 
appropriate. On both panels, the dotted vertical lines indicate the limits of the age bins. 



model • 
Siegler et al. (2007) o 
Corpenter et ol. (2009) □ 




Figure 7. Time evolution of the fraction of stars with 24^m 
excess for the best-fit model (black solid line and filled circles), 
compared with the findings of Siegler et al. (2007, blue diamonds) 
and Carpenter et al. (2009, red squares). The plot indicates that 
our model evolves on a longer timescale as found by these studies, 
although it is in agreement with the results of Siegler et al. (2007) 
for young (< 100 Myr) systems and Carpenter et al. (2009) for 
older ones. Note that for these studies, an excess was defined as 
R24 > 1.15, so we adopted this threshold for this plot. 



have 2 separate disk components: a component with r <3.5 
AU (rec ent work als o shows that for this component r > 0.5 
AU, e.g. Smith et al.. 2009. ) and a cold component at r-^^lSO 
AU (|Smith et all 1200^ ). with the middle region possibly 
cleared by a planetary system ijWvatt et al.ll2005h . Our re- 
sults suggest that the hotter disk component is transient, 
whereas the cold component is evolving in st eady state. 
HP 6 9830 has 3 known planetary companions (|Lovis et al.l 
l2006l). as well as an asteroid belt around 1 AU from the star 
I Beichman et al ] |2005l : iLisse et al.|[2007l ). Since a disk with 
such a small radius is expected to have processed its mate- 
rial at 4.7 Gyr (the age of HD 69830), this system can be 



considered anomalous; fitting it in our models would require 
unrealistic parameter values. 

HD 10647 also has a known Jupiter-mass planetary 
companion and a large cold disk at a radius of 300 AU has 
been detected at sub-mill imetre wavelengths (jLiseau et al.l 
l2008l ). lLawler et all |200i) fit the IRS spectrum with a disk 
radius of 16-29 AU, while |l iseau et all (|2010t ) imaged this 
disk with Herschel, finding a peak in the dust emission at 
85 AU. This system has the largest 70/.im excess in our sam- 
ple. As for HD 69830, finding a model that fits this object 
would require extreme parameter values, as it is unusually 
bright compared to other disks of similar age and colour 
temperature. This could be caused by the disk around HD 
10647 having unusually strong planetesimals, or having been 
recently stirred. 



5.2.2 Model predictions of24fim excess detection rates 

We compared the detection rates of excess emission at 24/im 
predicted by our models to each of the three surveys used in 
this analysis. Since at 70/im the large scatter of excesses and 
accuracy of flux measurements makes determining an excess 
threshold difficult and excess is quantified instead with the 
X70 statistic (Eg. I12|l . we only "predict" detection rates at 
24^m. We use the relevant age, spectral type and distance 
distribution of each survey to estimate their detection rates 
from the gen e rated model population. For the survey by 
iTrilling et al.l (|2008l ). we find a predicted detection rate of 
3.3%, in excellent agreement with their detection rate of 
4. 21^'?%. We also find a good agreement with the results 
of IBeichman et al] (|2006al 'l. who found a detection rate of 
7.3%, while our model p redicts 7.1%. Finally, the study of 
IHillenbrand et al] (|2008l ) constrained the 24/xm rate to be 
< 40%, with our model predicting a rate of 16.9%. 
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6 DISCUSSION 

6.1 Significance of best-fit parameters 

The best-fit values we find for Q}^ and Dc are significantly 
differen t from those found for the models of IWvatt et alj 
(|2007bl ) for debris disks around A stars. Both and Dc 
are an order of magnitude larger (values for these param- 
eters in the A stars study are 150 Jkg~^ and 60 km re- 
spectively). This could suggest either that our models are 
incomplete and these parameter values compensate for the 
value of one of several other free parameters being poorly 
chosen or fitted, or that the properties of debris disks around 
later-type stars are actually different. Other regions of pa- 
rameter space were also explored, and our best fit constitutes 
the best trade-off between the goodness-of-fit statistic for 
the evolution of excess rates (Fig. [S]and[S|), and realistic pa- 
rameter values, as stated in Sec. [S] 

In comparison with the study of A stars, we find dif- 
ferent values for rmin, 7, Afmid, and QJj. The most sig- 
nificant differences are the higher values of QJ, and Dc, al- 
though the exact values of the individual free parameters 
are poorly constrained and are not as informative as combi- 
nations of these parameters. In this case, the values point to 
a slower evolution of -F70, tot /-F70, phot and of F24,tot / F2i,phot 
compared to what was found for disks around earlier- type 
stars (e.g. lSu et al.|[200^ V It should also be emphasised that 
while FEPS data is available at several wavelengths for some 
of the objects presented here, we only consider data at 24 
and 70/^m. Using different wavelengths colours might result 
in different radii, temperatures and lumin osities, for exam- 
le if there is cool dust present in the disk. ICarpenter et al.l 
2OO9I ) indeed found that the dust temperatures derived 
from MIPS 24 and 70/im were lower than those inferred 
from spectra obtained with IRS (Infrared Spectrograph), and 
[Hillenbrand ct al,, (2008i ) find different colour temperatures 
when using 33 and 24/^m or 70 and 33/im fiux density ratios. 

Q*0 might vary because the composition of the disk may 
be different around later-type stars, or the bodies making 
up the disk may be more compacted by previous collisions, 
making them stronger and accounting for a higher intrinsic 
strength of the disks. The value of Dc indicates an initial 
population made of large objects, up to Pluto-size asteroids 
rather than smaller planetesimals, hinting at Kuiper-belt- 
like properties. However, the difference between solar-type 
and A stars results may also be due to inaccuracies in deriv- 
ing the radius of disks with the 24 — 70/im colour tempera- 
ture. 

Because small grains are inefficient emitters, and we are 
assuming blackbody radiation from all grains in the disk, the 
actual radii may be larger than the derived values. This is 
confirmed by results of disk imaging which show that disks 
are generally 2-3 ti mes larger than the blackbody value (e.g. 
[m aness et al.ll2009l '). We therefore emphasise that in this pa- 
per we are finding an effective value for the parameters 
and Dc, as was done for the A stars study. The discrepancy 
in effective values between A and FGK stars could then be 
caused by an dearth of small grains around A stars com- 
pared to solar-type stars. This is because the blowout size 
of grains is larger around earlier-type stars (~ 10 /xm for A 
stars and ~ 1/im for a GO star), meaning that with larger 
number of inefficient emitters, the emission is less well de- 
scribed by blackbody emission for disks around FGK stars. 



Hence the ratio of real radius to blackbody radius could be 
larger for solar-type stars than for A type stars. Despite ef- 
fective values being different, the real values of the planetes- 
imal strength and largest planetesimal size may therefore be 
similar for solar-type and A stars. In the next section, we 
discuss and quantify this. 



6.2 Effective and real parameters 

iBonsor fc WvattI (|2010l ) used grain modelling, taking into 
account the non-blackbody nature of the grains and a real- 
istic size distribution, to show that for A stars, real disk radii 
are expected to be larger than blackbody radii by a factor of 
~ 2 (with some dependence on radius and spectral type), in 
agreement with observations of resolved disks. They also de- 
rived scaling laws between realistic and effective values for 
parameters within the context of the modelling presented 
here. Assuming that the real radii are larger by a factor Xr 
compared to blackbody radii, then the real parameters can 
be calculated from parameters which were found using the 
blackbody assumption, by making sure that all disks have 
the same luminosity evolution when the radii are changed. 
In the following discussion we use a simplified notation so 
that M = Mmid, Q = Q*D and D = Dc- Furthermore, we 
use the subscript r to denote a real value, while its absence 
denotes an effective value: Qr refers to the real planetesimal 
strength, while Q refers to the fitted (effective) parameter. 
Disks start off with the same luminosity if 



MD" 



(13) 

from Eq. 3-4 of IWvatt et"aLl (|2007bh (see also Eq. 16 of that 
paper), and disks evolve on the same timescale if 



-)5/6 ^0.5 ^T/'.i _ 



q5/6 ^0.5 



(14) 



Note that these constraints are for an average popula- 
tion, and that they can be violated on a case-by-case basis. 
The real parameters for each population cannot be derived 
from the parameters of the blackbody fit, because there are 
4 free parameters and 2 constraints. However, it is possible 
to compare the real parameters of the FGK and A star pop- 
ulations from their blackbody fits. We rearrange Eq. (|13[l 
and (|14|) to derive expressions for Xr and Qr for the FGK 
stars population in terms of the equivalent A stars parame- 
ters. We use the additional subscripts A for A stars and F 
for FGK stars in the following comparison. 

If we assume that the strength of the planetesimals fol- 
lows Qr oc Dr'^ (i.e. that the strength of the planetesimals is 
mainly gravitational), and that planetesimals have the same 
composition around A and FGK stars, we can also rearrange 
these equations to fi nd that, with t he par ameter values found 
in this paper and in lWvatt eral] (|2007bt ) (Qa = 150 Jkg"\ 
Dc — 60 km and M = 10 M^) as we l l as t he value of Xr 
found for A stars bv lBonsor fc WvattI (|2010l '). 



r,F 



Mr. 



3.4 



r,A 



Dr 



(15) 



and using this to derive an expression for the scaling factor 
yields 
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~ 4.. 



Dr.A 
Dr,F 



3/4 



(16) 



We can expect Xr to be larger for FGK stars than for A 
stars for the reason mentioned earlier that there are more in- 
efficient emitters around FGK stars. Unfortunately the over- 
lap between the sample of disks resolved with Spitzer and 
the data sample used in this paper is too small for this to 
be useful to constrain X r for FGK stars, alth ough imaged 
disks such as HD 181327 (|Schneider et al.ll2006l ') suggest that 
a value of ~ 3 is a goo d estimate for Xr,F- Using Xr,A ~ 2 
(IBonsor fc Wvattll2010l ). the value of Xr f is consistent with 
the estimate of ^ 3 t hat is expected from imaged disks (e.g. 
[Schneider et al ] |2006l ) if the planetesimals around FGK stars 
are larger than those around A stars by a factor of ~ 5. 
Putting this number into Eq. (|15p . this means that the me- 
dian disk mass around FGK stars is ~ 0.72 times that for 
disks around A star s. This is in agreement with the findings 
of e.g. iNattal l|2004 ) that earlier-type stars are expected to 
have more massive disks. 

Finally, the values found for 7, rmin and Tmax point to 
a slightly flatter disk distribution for FGK stars compared 
to their earlier-type counterparts. This, however, does not 
take into account the possible systematic difference in actual 
disk radius contained in the scaling factor Xr,F, which would 
result in more large disks around FGK stars to what was 
found for the A stars. 



6.3 Time evolution of the 24/im statistics 



As pointed out in Section [5T2l the temporal evolution we de- 
rive in the model presented here is signifi cantl y slower than 
that o bserved by ICarpenter et al.l (|2009l ) and ISiegler et al.l 
The need for the model to yield slower time evolution 
is most obvious on Fig. |S] and 1101 the latter of which shows 
that several disks in the sample have small radii despite be- 
ing relatively old and therefore expected to have processed 
all their mass. 

We suggest two possible reasons for this discrepancy. 
In the first , the faster evoluti on seen by ICarpenter et al.l 
l|2009l ) and Isiegler et all (120071 ) could be due to the disks 
having two components. Surveys at 24/xm that focus on 
young stars would trace the hotter component that evolves 
on short timescales (being at smaller radii), whereas sur- 
veys of older stars would trace the colder component. In- 
deed, the data we use in this paper comes from surveys 
looking at older stars, for which the hot component would 
have faded below detec table levels. In contrast, the sample of 
ICarpenter et al. (l2009l) has o ver 80% of stars younger than 1 
Gvr. and lSiegler et al.l (|2007|) h a ve ove r 70% of young stars. 

Although Carpenter et al.l (|2009l l conclude from IRS 
spectra that 24/im excess emission arises from cold dust, 
the spectrum of a two-component disk could still look like 
a cold dust emission spectrum if the hot component, while 
bright enough to produce detectable 24^m excess emission 
at young ages, is too faint to make the 24 — 70/xm colour look 
like hot dust emission. Therefore more evidence is needed to 
determine whether two-component disks are responsible for 
the discrepancy in observed evolutionary timescales. 

The alt ernative explanation is that the fast 24/jm evolu- 
tion seen bv lCarpenter et all (|2009l ) and lSiegler et all (|2007l ') 



is the true evolution of single temperature disks (e.g. Kuiper 
belts), and that this really does disappear by a few 100 Myr 
(see Fig. I6I7|I . The 24/im emission from the older systems 
that contribute to the survey statistics modelled in this pa- 
per then must arise from a transient component. This is 
already thought to be the case for some of the old systems 
for which the 24pm e mission comes from a hot component 
within a few AU (e.g. IWvatt et aLll2007bl '). and is perhaps 
favoured by the fact that IRS spectra do not show evidence 
for two-component disks, but this in not the case for all of 
the older systems, such as HD 10647 which would have to be 
a significant outlier in such a model. 

Here we have modelled one of these possible cases, 
where 24/im comes from the cold component of a two- 
component disk. It is difficult to favour an interpretation 
with the small number of disks detected at both 24 and 
70^m, but additional detections from future surveys will 
help distinguish between these possible explanations. Fi- 
nally we also note that as in this paper we use a single- 
temperature disk model; fitting this data with an extended 
disk model and a range of temperatures could also signifi- 
cantly improve the fit. 



6.4 Fractional luminosity vs. radius 

Fig. [S] shows a plot of fractional luminosity / against disk 
radius, with the observed data plotted over the model pop- 
ulation. Also plotted are the detection limits at both 24 and 
70/im, given by Eq. (|10p . and the lines of maximum frac- 
tional luminosity at 2 and 10 Gyr, given by Eq. for 
an FO dwarf and the best-fit model parameters. The detec- 
tion thresholds are calculated assuming calibration limits of 
R24 ~ 0.11 and R70 = 0.5. In theory, disks should lie on 
or below the line of maximum lumin osity for their age an d 
spectral type, although as noted by IWvatt et al.l (|2007bD . 
the precise location of these lines depends on parameters 
which may vary between disks such as QJj, e, Dc, and the 
spectral type of the disk's host star. The sharp increase in 
the 24/im threshold for disks with radii larger than ~ 50 
AU means that the disks with larger radii, i.e. those disks 
which are expected to have the strongest emission at 70/im 
will be more difficult to detect at both wavelengths as they 
will require high fractional dust luminosity. 

Values of ///max for disks detected at 24 and 70^m are 
given in Table (2] Since /max depends on the disk radiu^, a 
good estimate can only be made for those disks detected at 
both wavelengths, so we do not include a value ///max for 
disks detected at 70/im only. Several disks st and out in Ta- 
ble [2| and on Fig.[8l In particular, HD 69830 (jTrilling et al.1 
[2003) has a luminosity / over 800 times the maximum theo- 
retical value /max for its age, radius and spectral type. This 
comes from a very small disk radius of 1 AU combined with 
a relatively old stellar age of 4.7 Gyr. A fit to the model pop- 
ulation yields / oc r"'^^, while a fit to the observed popula- 
tion of disks in the lTrilling" et al.l (|2008h and lBeichman et al.l 



^ Note that the calculation of /max is unaffected by the discussion 
in Sec. l6.2l on the difference between real and effective radius, as 
long as /max is calculated using both effective parameters and 
blackbody radius (which would give the same value that would 
be calculated from the real parameters and real radius). 
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Figure 8. Dust luminosity / = Lij^/Lt plotted as a function of disk radius -Rdust ■ The left-hand panel shows the whole model population, 
while the right-hand panel only has the model population that could be detected at 24 and 70/^m. The model population is shown with 
small dots, and data symbols are the same as for Fig. [6] Detection thresholds for an FO dwarf arc indicated by a solid (70^m) and a 
dashed (24/im) line. Lines showing the maximum possible fractional luminosity for an FO dwarf are shown as dashed lines for 2 and 10 
Gyr. The dotted line indicates the limit, given by Eq. Hilt , where Poynting-Robertson drag becomes important for an FO dwarf. 



(|2006al ) samples, without anomalously bright systems, and 
without disks known to have two components (HD 109085), 
gives a relation / oc ;-0-6S±o.05^ excellent agreement. 



6.5 Fractional luminosity vs. age 

Fig. [9] shows the fractional luminosity of disks as a function 
of age, as well as theoretical evolution lines for disks of initial 
mass 3, 30 and 300 of radii 10 and 100 AU. The larger 
disks do not reach steady state in the timescale considered 
here, while the smaller disks have reached their coUisional 
equilibrium by 1 Gyr. 

A fit to the model population gives the relation / oc 
f^'^'*^ while a fit to the observed population yields / oc 
^-o.i8±o.io^ just within 3 a of the best-fit model. Fitting 
of the observed disks is strongly affected by the low num- 
ber of disks in our data set, especially by the relatively few 
young disks (< 1 Gyr) in the sample, but one can conclude 
from Fig. [9] that the model population is a reasonable fit 
for the observed population. The validity of these fits will 
benefit from additional observations, which will allow more 
significant conclusions to be drawn with larger samples. In- 
teresting to notice is the absence of a peak in excess flux 
analogous to the o ne seen in the obs erved A star samples 
around 10-15 Myr (|Currie et al.ll2008l '). 



6.6 Radius vs. age 

Fig.[TD]is a plot of disk radius as a function of age, for those 
disks that could be detected at 24 and 70/im. 

A lower limit for detectable radii is visible in the model 
population, with this lower limit increasing with age beyond 
~100 Myr. Older disks must therefore have larger radii to 
be detected at both 24 and 70/im. This comes from the fact 
that disks with larger radii are slower in processing their 
mass and can therefore remain above detection thresholds 
for a longer time than disks with small radii. This radius in- 
crease with age is therefore only due to detection sensitivity 
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Figure 10. Disk radius plotted as a function of age. The model 
population is shown with small dots, and data symbols are the 
same as for Fig. \E\ Only the populations which could be detected 
at both 24 and 70 /im are plotted. Also shown are power-law fits 
to the model (dashed line) and observed (solid line) populations. 
These are not in agreement, and possible reasons for this are 
discussed in the text. 



considerations. In fact, there is a s yet no evidence fo r a cor- 
relation between radius and age (jNaiita fc Williamsll2009' ). 

A fit to the model population confirms the apparent in- 
crease of disk radius with age, with the relation r oc t^'^^, 
while a fit to the observed population yields r oc ^o ofiO M, 
The fit to the d ata is str o ngly infiuenced by the disks 
from the [Hillenbrand et al.l (|2008l ) sample, which consists 
of younger, and mostly larger disks. We also see a lack of 
observed small and young disks, which are predicted by the 
model. Although this could be evidence for rapid inner clear- 
ing, the small number of disks of the observed sample makes 
it difficult to draw robust conclusions as to the cause of this. 
Thus we note that it is worth scrutinising <300 Myr stars to 
see if there really is an absence of <10 AU disks that would 
be the precursors of older systems such as HD 40136 and HD 
109085 {ri Corvi), and whether the radii derived from 24-70 
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Figure 9. Dust luminosity / = Ljr/L, plotted as a function of age. The model population is shown with small dots, and data symbols 
are the same as for Fig. [B] For < 3cr detections, the fluxes plotted are maximum values. Left: Entire model population and observed 
objects detected at both 24 and 70/im. Also shown are theoretical evolutionary tracks for disks of radii 10 (solid line) and 100 (dashed 
line) AU, and for initial disk mass of 3, 30 and 300 (with the higher masses corresponding to the lines with higher early-age fractional 
luminosities), around an FO star. Right: Same with only the model population that could be detected at 24 and 70fim plotted. Also 
plotted on the right-hand panel are power-law fits to the model (dashed line) and the data (solid line). The power law fit to the model 
population is in agreement with the fit to the data (see text). 



colour temperatures appear smaller than reality for a 
significant number older disks, perhaps due to the presence 
of both ho t and cold dust (as is known to be the case for 
Corvi; see IWvatt et al1l2005l l. or due to the action of P-R 
drag (dis cussed later). In comp arison, for sources detected 
at 70/im, ICarpenter et al.l (|2009l ) fitted FEPS IRS spectra to 
derive probability distributions for inner disk radii peaking 
at a few tens of AU. 

The obs e rved disks HD 69830 and HD 166 
iTrilling et al.l |2008| ) appear on the plot where the 
model predicts that no disks should be detected (bottom 
right corner of Fig. \W\ For HD 69830, we attribute this to 
an unusually high disk luminosity; indeed this disk is found 
to have over 800 times the maximum theoretical luminosity 
for its radius, age and the spectral type of their host, as 
was already mentioned in Sec. 16.41 HD 166 lies just outside 
the model population. Its disk is small for its age and its 
luminosity is just over its ma^ximum theoretical luminosity; 
the unusual brightness for a disk of this age and colour 
temperature may mean that the disk around HD 166 has 
unusually strong or massive planetesimals compared to its 
peers. Other disks are also, although to a lesser extent, 
anomalously bright for their age: HD 109085, HD 10647 and 
HD 207129. Unusual disk properties might explain this; HD 
109085 has two disk components, as discussed above, while 
recen t imaging of the disk around HD 207129 (|Krist et al.l 
l2010h and HD 10647 (Stapelfeldt et al., in preparation) 
with the ACS coronagraph on the Hubble Space Telescope 
shows that these disks have radii larger than the 24-70/im 
radii used here by an order of magnitude. We discussed in 
Sec. I6.2l whv we might expect blackbody radii to be smaller 
than the true radii, and these discrepancies also indicate 
that the 24-70/xm radius could be a worse estimator for 
the true radius for FGK stars, compared to A stars. This 
could be because some of the disks have two components, 
as discussed in Sec. 16.31 or because there is a large spread in 
the Xr,F factors between the disks. Larger disk samples and 



detailed analysis should help improve our understanding of 
this efi'ect. 

Dust might also find itself closer to the star than the 
planetesimal belt under the action of Poynting-Robertson 
drag. This becomes important when P-R drag and collisional 
timescales are comparable, or in terms of luminosity, when 
/ < /pR, with /pR, being given by Eq. ([TT}. 

This limit is shown on Fig. [8] From this plot, it emerges 
that disks with small radii could be affected by PR drag, 
given the proximity of the limit for PR drag to the limit of 
detection at 24^m for these disks. In our sample, no sys- 
tem has ///pR < 1, meaning that PR drag is not significant 
for the observed disks we are considering ; the most likely 
system to be affected by this effect is HD 40136. We also 
note that if the disk radius has been underestimated by a 
factor Xr , then the limit at which P-R drag becomes impor- 
tant given in Eq. (|lip should be lower by a factor Xr^^^- 
Therefore we do not expect P-R drag to have significantly 
influenced the blackbody radius, derived from dust temper- 
ature, which is thus a good proxy for the location of the 
planetesimal belt (with the caveats mentioned in Sec. 16. 2p . 

6.7 Histogram of ///max 

Fig. 111! shows a histogram of the quantity ///max, plotted 
for our data sets and for the model. We also plot the model 
with an arbitrary Gaussian smoothing of 1 dex to account 
for variations of disk properties between individual systems, 
which we have not considered in our model. The range of 
fluxes i s broader in the data, which, as stated bv lWvatt et al.l 
(|2007bD can be explained by the fact that we have assumed 
in our model that all the planetesimal belts have the same 
properties, where in reality the values of parameters such as 
Dc, Qh and e are expected to vary from one belt to another. 
The parameters we used would then be correct as an average 
over a whole population of stars, but not necessarily for in- 
dividual debris disks. This is particularly relevant when the 
sample of observed disks is as small as it is here. Although 
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Figure 11. Histogram of dust fractional luminosity for the ob- 
served population with an excess at both 24 and 70 /im (solid 
line), the model population (dashed, red) and the model popu- 
lation with a Gaussian smoothing of 1 dex (dotted, blue). Lumi- 
nostiies are expressed as a fraction of the maximum theoretical 
luminosity /max , given by Eq. JgJ , and only the model population 
which could be detected at both 24 and 70 ^m is plotted. 



observed in these surveys are made up in large part of young 
(< 1 Gyr) stars, around which the hot component would still 
be bright enough to produce detectable amounts of 24/im 
emission. However, this excess emission would evolve rapidly 
due to the disk component's small distance from its host 
star. In contrast, the data we use in this analysis includes 
many older stars, for which the 24/xm emission would come 
from the cold component of the disk, which would evolve 
on longer timescales. However, an alternative explanation 
for the apparent slow evolution might also be that it comes 
from cold dust, with the excess emission past a few Myr then 
caused by transient events. 

Small data samples are strong caveats on these results. 
In most cases, the flux statistics we derive from our analysis 
are based on only a few disks. This shows the importance of 
collecting more observations of disks around late-type stars 
with large-scale surveys. More observations would help de- 
termine why the 24 — 70/im radii might be worse radius esti- 
mators around stars of different types, and whether the slow 
evolution of 24/im flux statistics we derived in this paper is 
statistically significant. 



our data set is small, main features are reproduced by the 
smoothed model. 



7 CONCLUSIONS 

In this paper, we have showed that the main features seen 
in observations of debris disks around FGK stars can be 
attributed to collisional grinding of planetesimals. We mod- 
elled data collected by teams with the Spitzer Space Tele- 
scope, and used 24 and 70/im colours to derive simple prop- 
erties of the disks for each of these systems. Our results are 
consistent with models of planetesimal strengths and their 
size distributions, although we find parameters correspond- 
ing to these properties to be an order of magnitude difi'erent 
from those found in the study of A st ars that was done using 
the same model (|Wvatt et al.ll2007bl ). 

We discussed whether these values might be effective 
values due to differences in stellar environment, and found 
that properties of disks around later-type stars may be some- 
what different from those around earlier-type stars, because 
they evolve on longer timescales and the refore have had more 
time to form larger planetesimals (e.g. iKennedv fc KenvonI 
[iooi), which are then also stronger if the strength is mainly 
due to gravitational pull. We also find that disks around 
FGK stars are less massive than those around A sta rs by a 
small factor, in agreement with previous findings (e.g. lNat"tal 
[iooi). Future observations to increase the sample of debris 
disks detected at multiple wavelengths will allow us to test 
our assumption that the blackbody disk radii underestimate 
the true radii of disks around FGK stars by a larger fac- 
tor than for those around A stars. Realisti c grain modelling 
analog ous to the work done for A stars bv lBonsor fc WvattI 
l|2010f ) could also help constrain this factor. 

We proposed two explanations for the discrepancy be- 
tween the slower time evoluti on we derive in this a naly- 
sis compared to th at found by ICarpenter et al. I (|2009l ) and 
ISiegler et all (HqO^). We first suggested that this could be 
due to disks having two components. The samples of disks 
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